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Ab initio theoretical calculations are reported for the electric (El) dipole allowed and inter- 
j> ! combination fine structure transitions in Fe V using the Breit-Pauli R-matrix (BPRM) method. 
' We obtain 3865 bound fine structure levels of Fe V and 1.46 x 10^ oscillator strengths, Ein- 

o ; 

^ ' stein A-coefficients and line strengths. In addition to the relativistic effects, the intermediate 

o ■ 

^ coupling calculations include extensive electron correlation effects that represent the complex con- 

I figuration interaction (CI). For bound-bound transitions the BPRM method, based on atomic 

^ collision theory, entails the computation of the CI wavefunctions of the atomic system as an 
>^ 

■ (electron + target ion) complex. The target ion Fe VI is represented by an eigenfunction expan- 
J> ', sion of 19 fine structure levels dominated by the spectroscopic configuration 3d^, and a number 
^ of correlation configurations. Fe V bound levels are obtained with angular and spin symmetries 
■ - - ' SLn and Jn of the (e + Fe VI) system such that 25"+ 1 = 5,3,1, L < 10, J < 8. The bound levels 
are obtained as solutions of the Breit-Pauli (e + ion) Hamiltonian for each Jn, and are designated 
according to the 'collision' channel quantum numbers. A major task has been the identification 
of these large number of bound fine structure levels in terms of standard spectroscopic desig- 
nations. A new scheme, based on the analysis of quantum defects and channel wavefunctions, 
has been developed. The identification scheme aims particularly to determine the completeness 
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of the results in terms of all possible bound levels with n < 10,1 < n — 1, for applications to 
analysis of experimental measurements and plasma modeling. Sample results are presented and 
the accuracy of the results is discussed. A comparison of the dipolc length and velocity oscillator 
strengths is presented, indicating an uncertainty of 10-20% for most transitions. 

1. Introduction 

Transition probabilities of heavy elements, particularly the iron group, are of great importance 
in astrophysical and laboratory sources. Fuhr et al. [1] have compiled data from a number of 
available sources. However, the accuracy and the extent of these data is largely inadequate for 
many general applications such as the calculation of local thermodynamic equilibrium (LTE) 
stellar opacities [2,3], and radiative levitation and accelerations of heavy elements [4]. Among 
the particular applications including Fe V as a prominent spectral constituent are the the non- 
LTE models of Fe V spectra in hot stars [5], and the observed extreme ultraviolet Fe V emission 
from young white dwarfs [6]. For example, currently available data for Fe V fails to account 
for the observed opacity of iron in the XUV region where observations of newly formed hot and 
young white dwarfs clearly show Fe V hues [6] . In all of these apphcations it is highly desirable 
to have as complete a dataset of radiative transition probabilities as possible. While the twin 
problems of completeness and accuracy pose a challenge to the theoretical methods, they are of 
interest not only in various applications but may also be of use in the analysis of experimental 
measurements of observed energy levels of complex atomic systems from the iron group. 

The Opacity Project (OP) [7,2] and the Iron Project (IP) [8] laid the foundation for large- 
scale theoretical calculations using ab intio methods. The R-matrix method [9], based on atomic 
collision theory techniques and adapted for the OP [10] and the IP [8], has proven to be very 
efficient for these calculations. Whereas the OP calculations were all in the LS coupling approx- 
imation, with no relativistic effects included, the subsequent IP work is in intermediate coupling 
using the Breit-Pauli extension of the R-matrix method [8]. While most the IP work has con- 
centrated on collisional calculations, recent works have extended the BPRM method to radiative 
bound-bound and bound-free calculations for transition probabilities [11], photoionization [12], 
and (electro- ion) recombination [13]. The first comprehensive BPRM calculation of fine structure 
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transition probabilities was carried out for the highly charged ions Fe XXIV and Fe XXV [11] 
that are of special interest in X-ray astronomy. Very good agreement was found with existing 
results available for a limited number of transitions but using very accurate theoretical methods 
including relativistic and QED effects [14,15], thus establishing the achievable accuracy for the 
BPRM calculations. However those He-like and Li-like atomic systems are relatively simple, and 
the electron correlation effects relatively weak, compared to the low ionization stages of iron 
group elements. The present work attempts to enlarge the scope of the possible BPRM calcula- 
tions to include the iron group elements, as well as to solve some outstanding problems related 
to level identifications in ab initio theoretical calculations using collision theory methods. 

Unlike atomic structure calculations, where the electronic configurations are pre-specified 
and the levels identified, the bound levels calculated by collision theory methods adopted in the 
OP and the IP need to be identified since only the channel quantum numbers are known for 
the bound states corresponding to the (e -|- ion) Hamiltonian of a given total angular and spin 
symmetry SLn or Jtt. The precise correspondence between the channels of the collision complex, 
and the bound levels, must therefore be determined. The problem is non-trivial for complex 
atoms and ions with many highly mixed levels due to configuration interaction. In the OP work, 
carried out in LS coupling, this problem was solved by an analysis based on quantum defects and 
the numerical components of wavefunctions in the region outside the R-matrix boundary (that 
envelops the target ion orbitals) . The present work extends that treatment to the analysis of fine 
structure levels computed in intermediate coupling. In addition, considerable effort is devoted 
to the determination of the completeness of the set of computed bound levels; comparison with 
the expected levels derived from all possible combination of angular and spin quantum numbers 
reveals the missing levels. The general procedure could be applied to spectroscopic measurements 
and the analysis of observed levels of a given atomic system by comparison with the theoretical 
predictions. 

2. Theory 

The general theory for the calculation of bound states in the close coupling (CC) approximation 
of atomic collision theory, using the R-matrix method, is described by Burke and Seaton [16] and 
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Seaton [17]. The application to the Opacity Project work is described by Seaton [7], Berrington 
etal. [10], and Seaton etal. [2]. The relativistic extensions of the R- matrix method in the Breit- 
PauH approximation are discussed by Scott and Taylor [18], and the computational details by 
Berrington, Eissner, and Norrington [19]. The application to the Iron Project work is outlined 
in Hummer etal. [8]. 

In the present work we describe the salient features of the theory and computations as they 
pertain to large-scale BPRM calculations for complex atomic systems. Identification of fine 
structure energy levels is discussed in detail. 

Following standard collision theory nomenclature, we refer to the (e + ion) complex in terms 
of the 'target' ion, with N bound electrons, and a 'free' electron that may be either bound or 
continuum. The total energy of the system is either negative or positive; negative eigenvalues of 
the (N + l)-electron Hamiltonian correspond to bound states of the (e + ion) system. In the 
coupled channel or close coupling (CC) approximation the wavefunction expansion, "^{E), for 
a total spin and angular symmetry SLn or Jn, of the (N+l) electron system is represented in 
terms of the target ion states as: 

^^Eie + ion) = AY,Xi{ion)ei + Y,Cj<^j, (1) 

i 3 

where Xi is the target ion wave function in a specific state SiLiiTi or level JjTTj, and 9i is the wave 
function for the (N-l-l)th electron in a channel labeled as SiLi(Ji)7ri kf£i{SLn) [Jtt]; kf{— e^) is 
the incident kinetic energy. In the second sum the <l>j's are correlation wavefunctions of the (N+l) 
electron system that (a) compensate for the orthogonality conditions between the continuum and 
the bound orbitals, and (b) represent additional short-range correlation that is often of crucial 
importance in scattering and radiative CC calculations for each SLn. 

The functions ^(E) are given by the R-matrix method in an inner region r < a. These are 
bounded at the origin and contain radial functions that satisfy a logarithmic boundary condition 
at r = a [20]. In the outer region r > a the inner region functions are matched to a set of linearly 
independent functions that correspond to all possible (e -|- ion) channels of a given symmetry 
SLn or Jn. The outer region wavefunctions are computed for all channels, {CtStLtnt)el, where 
Ct is the target configuration, and used to determine the individual channel contributions (called 
"channel weights"). 
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In the relativistic BPRM calculations the set of SLtt are recoupled to obtain (e + ion) levels 
with total Jtt, followed by diagonalisation of the (N+l)-electron Hamiltonian, 

H^l^^ = E^. (2) 

The BP Hamiltonian is 

rj-BP TT I rrmass i rrDar , rrso (q\ 

where i?iv+i is the nonrelativistic Hamiltonian, 

N+l ( 9 7 N+1 r, 1 

i/.«=E -v?-2^+i:f . (4) 

and the additional terms are the one-body terms, the mass correction term, the Darwin term 
and the spin-orbit term respectively. Spin-orbit interaction, splits the LS terms into fine- 

structure levels labeled by Jtt, where J is the total angular momentum. Other terms of the 
Breit-interaction [22] , 

H"" = Ebi,- (so + so') + (ss')] , (5) 

representing the two-body spin-spin and the spin-other-orbit interactions are not included. 

The positive and negative energy states (Eq. 1) define continuum or bound (e -|- ion) states, 

E = > — >• continuum (scattering) channel , . 

E = < — ^ bound state, 

where p is the effective quantum number relative to the core level. If £" < then all continuum 

channels are 'closed' and the solutions represent bound states. Determination of the quantum 

defect i/J'ii)), defined as i^j = n — //(£) where Uj, is relative to the core level SiLiJTi, is helpful in 

estabhshing the £- value associated with a given channel (level) . 

At E < a scattering channel may represent a bound state at the proper eigenvalue of the 

Hamiltonian (Eq. 2). A large number of channels arc considered for the radiative processes of 

Fe V. Each SLtt or Jtt symmetry is treated independently and corresponds to a large number 

of channels. Therefore, the overall configuration interaction included in the total (e -|- ion) 

wavefunction expansion is quite extensive. This is the main advantage of the CC method in 

representing electron correlation accurately. 



a) Level identification and coupling schemes 

The BPRM calculations in intermediate coupling employ the pair-coupling representation 

Si ~\~ L-i ^ J -i 

J, + i^K , (7) 
K + s — > J 

where the 'i' refers to the target ion level and £, s are the orbital angular momemtum (partial 
wave) and spin of the additional electron. According to designations of a collision complex, a 
channel is fully specified by the quantum numbers 

{SiLi Ji)Tii eiUK s [Jtt] (8) 

The main problem in identification of the fine structure levels stems from the fact that the 
bound levels are initially given only as eigenvalues of the (e + ion) Hamiltonian of a given sym- 
metry Jtt. Each level therefore needs to be associated with the quantum numbers characterizing 
a given collision channel. Subsequently, three main parameters are to be determined: (i) the 
parent or the target ion level, (ii) the orbital, effective and principal quantum numbers (/, n) 
of the (N-l-l)th electron, and (iii) the symmetry, SLt:. The task is relatively straightforward for 
simple few-electron atomic systems. For example, in a recent work Nahar and Pradhan [11] have 
calculated a large number of transition probabilities for Li-hke Fe XXIV and He-like Fe XXV, 
where the problem of level identification is trivial, compared to the present work, since the bound 
levels are well separated in energy and in v. However when a number of mixed bound levels fall 
within a given interval {v, -|- 1), for the same Jtt, the quantum numbers and the magnitude 
of the components in all associated channels must be analysed. A scheme for identification of 
levels is developed (discussed later) that rests mainly on an analysis of quantum defects of the 
bound levels and their orbital angular momenta, and the percentage of the total wavefunction 
in all channels of a given Jtt. 

Following level identification, further work is needed to enable a direct correspondence with 
standard spectroscopic designations that follow different coupling schemes, such as between LS 
and JJ, appropriate for atomic structure calculations as, for example, in the NIST tables of 
observed energy levels [1] . The correspondence provides the check for completeness of calculated 
set of levels or the levels missing. The level identification procedure involves considerable ma- 
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nipulation of the bound level data and, although it has been encoded for general applications, 
still requires analysis and interpretation of problem cases of highly mixed levels that are difficult 
to identify 

b) Oscillator strengths and transition probabilities 

The oscillator strength (or photoionization cross section) is proportional to the generalized line 
strength defined, in either length form or velocity form, by the equations 

2 

(9) 
and 

2 

(10) 

In these equations lo is the incident photon energy in Rydberg units, and "Jfi and !P/ are the wave 
functions representing the initial and final states respectively. The boundary conditions satisfied 
by a bound state with negative energy correspond to exponentially decaying partial waves in all 
'closed' channels, whilst those satisfied by a free or continuum state correspond to a plane wave 
in the direction of the ejected electron momentum k_ and ingoing waves in all open channels. 

Using the energy difference, Ej^^ between the initial and final states, the oscillator strength, 
/ij-, for the transition can be obtained from S as 

/« - I^S, (11) 

and the Einstein's A-coefficient, Aji, as 

Ma.u.) = \a'^E%U,, (12) 

where a is the fine structure constant, and gi, gj are the statistical weight factors of the initial 
and final states, respectively. In terms of c.g.s. unit of time, 

Aj,is-') = (13) 

To 

where Tq — 2.4191~^^s is the atomic unit of time. 
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3. Computations 



The target wavefunctions of Fe VI were obtained by Chen and Pradhan [21] from an atomic 
structure calculation using the Breit-Pauli version of the SUPERSTRUCTURE program [22], 
intended for electron collision calculations with Fe VI using the Breit-Pauli R-matrix method. 
Present work employs their optimized target of 19 fine structure levels [21] corresponding to 
the 8-term LS basis set of 3d^(^F, ^P, ^G, ^p, 2^2, '^H, ^F, ^Dl). The set of correlation 
configurations used were Ss'^Sp^ScPAs, 3s'^3p^3(PAd, 3s3p^3d\ 3p^3d^, 3s'^3p'^3d^, and 3p^3dHs. 
The values of the scahng parameter in the Thomas-Fermi potential for each orbital of the target 
ion are given in Ref. [21]. Table I lists the 19 fine structure energy levels of Fe VI used in the 
eigenfunction expansion where the energies are the observed ones. Most bound levels in low 
ionization stages correspond to the level of excitation of the parent ion involving the first few 
excited states. The criterion remains the accuracy of the target represetation that constitute the 
core ion states. The (N-l-1) electron configurations, which meet the orthogonality condition 
for the CC expansion (the second term of the wavefunction, Eq. (1)) arc given below Table 1. 
The same set of configurations is used for all the states considered in this work. STGl of the 
BPRM codes computes the one- and two-electron radial integrals using the one-electron target 
orbitals generated by SUPERSTRUCTURE. The number of continuum basis functions is 12. 

The present calculations are concerned with all possible bound levels with n < 10, £ < n — 1. 
These correspond to total (e + Fe VI) symmetries {SLn) with (25" + 1) = 1,3,5 and L = - 10 
(even and odd parities). The intermediate coupling calculations are carried out on recoupling 
these LS symmetries in a pair-couphng representation, Eq. 6, in stage RECUPD. The computer 
memory requirement for this stage is the maximum, since it carries out angular algebra of 
dipole matrix elements of a large number of fine structure levels. The (e -|- Fe VI) Hamiltonian 
is diagonalized for each resulting Jtt in STGH. The negative eigenvalues of the Hamiltonian 
correspond to the bound levels of Fe V, that are found according to the procedure described 
below. 
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a) Calculation of bound levels 

The eigenenergies of the Hamiltonian for each Jtt are determined with a numerical search on an 
effective quantum number mesh, with an interval Azv, using the code STGB. In the relativistic 
case, the number of Rydberg series of levels increases considerably from those in LS coupling due 
to sphtting of the target states into their fine structure components. This results in a large number 
of fine structure levels in comparatively narrow energy bands. A mesh with A(z/) = 0.01 is usually 
adequate to scan for LS term energies; however, it is found to be of insufficient resolution for fine 
structure energy levels. The mesh needs to be finer by an order of magnitude, i.e., A(z^) = 0.001, 
so as not to miss out any significant number of bound levels. This considerably increases the 
computational requirements for the intermediate coupling calculations of bound levels over the 
LS coupling case by orders of magnitude. The calculations take up to several CPU hours per 
Jtt in order to determine the corresponding eigenvalues. All bound levels of total J < 8, of both 
parities, are considered. However, a further search with an even finer /S.v reveals that a few levels 
are still missing for some Jtt symmetries. 

b) Procedure for level identification 

The energy levels in the BPRM approximation (from STGB) are identified by Jtt alone. This 
is obviously insufficient information to identify all associated quantum numbers of a level from 
among a large set levels for each Jtt, typically a few hundred for Fe V. A sample set of energy 
levels for J = 2, even parity, obtained from the BPRM calculations is presented in Table II. 
The table shows energies and effective quantum number Vg, as calculated relative to the ground 
level (3(i^ ^i^3/2) of the core ion Fe VI. The complexity of the calculations, and that of level 
identification, may be gauged from the fact that 30 of these levels have nearly the same Ug. 
Further, the Vg do not in general correspond to the actual effective quantum number of the Fe V 
level since it may belong to an excited parent level, and not the ground level, of Fe VI. 

A scheme has been developed to identify the levels with complete spectroscopic information 
consisting of 

{Ct St Lt Jt Tidm J^, (14) 
and also to designate the levels with a possible SLn symmetry. The designation of the SLn, 
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from the identifications denoted above, is generally ambiguous since the collision channels are all 
in intermediate coupling. However, in most cases we are able to carry through the identification 
procedure to the LS term designation. An advantage of identification is that it greatly facilitates 
the completeness check for all possible LS terms and locate any missing levels. A computer code 
PRCBPID has been developed to identify all the quantum numbers relevant to the Jtt and the 
LS term assignments. Identification is carried out for all the levels belonging to a Jtt symmetry 
at a time. 

The components of the total wavefunction of a given fine structure energy level span all closed 
"collision" channels {CtStLt{Jt)'n't)el. Each channel contains the information of the relevant core 
and the outer electron angular momentum. The "channel weights" , mentioned earlier, determine 
the magnitude of the wavefunction in the outer R-matrix region of each channel evaluated in 
STGB. A bound level may be readily assigned to the quantum numbers of a given channel 
provided the corresponding channel weight (in percentage terms) dominates the other channels. 

The number of channels can be large especially for complex ions. For Fe V, for example, each 
level with J > 2 corresponds to several hundred channels. As the first step in the level identi- 
fication scheme we isolate the two most dominant channels by comparing all channel percentage 
weights. The reason is that the largest channel percentage weight may not uniquely determine 
the identifications since the channel weights are evaluated from the outer region contributions 
(r > a); the inner region contributions are unknown. Also, many levels are often heavily mixed 
and no assignment for the dominant channel may be made. 

The program, PRCBPID, sorts out the duphcate identifications in all the levels of the Jtt 
symmetry. Two levels with the same configuration and set of quantum numbers can actually 
be two independent levels due to outer electron spin addition/ subtraction to/from the parent 
spin angular momentum, i.e. St ± s — S. The identical pair of levels are tagged with positive 
and negative signs indicating higher and lower multiplicity respectively. The lower energies are 
normally assigned with the higher spin multiplicity. However, the energies and effective quantum 
numbers {i/) of levels of higher and lower spin multiplicity can be very close to each other, in 
which case the spin multiplicity assignment may be uncertain. 

One important identification criterion is the analysis of the quantum defect, fx, or the effective 
quantum number, v, of the outer or the valence electron. The principle quantum number, n, 
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of the outer electron of a level is determined from its u, and a Rydberg series of levels can be 
identified from the effective quantum number. Hence, in the identification procedure, u of the 
lowest member (level with the lowest principal quantum number of the valence electron) of a 
Rydberg series is determined from quantum defect analysis of all the computed levels for each 
partial wave I. The lowest partial wave has the highest quantum defect. A check is maintained to 
differentiate the quantum defect of a 's' electron with that of an equivalent electron state which 
has typically a large value in the close coupling calculations. The principle quantum number, 
n, of the lowest member of the series is determined from the orbital angular momentum of the 
outer electron and the target or core configuration. Once u and /i — n — i/ oi the lowest member 
are known, the n-values of all levels can be assigned for each paritial wave, I. The relevant 
Rydberg series of levels is also identified from the levels that have the same symmetry, Jtt, core 
configuration, Ct St LtTTt and outer electron orbital angular momentum /, but different u that 
differs between successive levels by ~ 1. While the ^{n £) are more accurate for the higher 
members of the series, they are more approximate for the lowest ones. The quantum defect of a 
given partial wave £ also varies slightly with different parent core levels and final SLJ symmetries. 

Of the two most dominant channels the proper one for each bound level is determined based 
on several criteria. There are cases when more than two levels are found to have identical 
identifications. These levels arc checked individually for proper identification. Often a swap of 
identifications is needed between the two sets of dominant channels since the second dominating 
channel is more likely to be associated with the given level, consistent with all other criteria. 
In some cases the most dominant channel (largest percentage weight in the outer region) may 
correspond to comparatively larger u for the partial wave i, than to a reasonable v for the second 
channel, indicating that the identification should correspond to the second channel. 

In a few cases a level is found not to correspond to any of the two dominant channels, 
predetermined from the channel weights. At the same time often a level is found to be missing in 
the same energy range. In such case the level is assigned to a channel of lower percentage weight 
that has a reasonable core configuration and term, nl quantum numbers for the outer electron 
and effective quantum number that match the missing level. 

There are a number of levels belonging to equivalent-electron configurations and require 
different identification criteria from those of the Rydberg states. These levels usually have: (i) a 
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number of approximately equal channel weights, and (ii) quantum defects that are larger than 
that of the lowest partial wave, or an inconsistent v that does not match with any reasonable 
n. Once these levels arc singled out, they are identified with the possible configurations of the 
core level, augmented by one electron in the existing orbital sub-shell. These low-lying levels are 
often assigned to those identified from the small experimentally available set of observed levels. 
The levels that can not be identified in the above procedure, such as by swapping of channels, 
or maching to a missing level, are assumed to belong to mixed states. These are not analysed 
futher by quantum defects. 

Two additional (and related) problems, as mentioned above, are addressed in the identifica- 
tion work: (A) standard LS coupling designation, S'Ltt, and (B) the completeness check for the 
set of all fine structure components within an LS multiplet. Identification according to collision 
channel quantum numbers is not quite sufficient to establish a direct correspondence with the 
standard spectroscopic notation employed in atomic structure calculations, or in the compiled 
databases such as those by the U.S. National Institute for Standards and Technology (NIST). 

The possible set of S'Ltts of a level is obtained from the target term, StLti^t, and the valence 
electron angular momentum, Z, at the first occurance of the level in the set. The total spin 
multiplicity of the level is defined according to the energy level position as discussed above. For 
example, the core 3(i^(^F'^) combining with a 4d electron forms the terms ^(P, D, F, G, HY and 
^{P^ F,G, Hy (Table IV) where the quintet for each L should be lower than the triplet. 

To each LS symmetry, SLn, of the set belongs a set of predetermined J-levels. The set of 
total J-values of same spin multiplicity is then calculated from all possible LS terms, equal to 
|L + S|. The program sorts out all calculated fine stucture levels with the same configuration, 
but with different sets of Jj and J, e.g. (C^ St Lf Jt TTtui) J % (including the sign for the upper or 
lower spin multiplicity), compares them with the predetermined set, and groups them together. 
Thus a correspondence is made between the set of SLn and the calculated fine structure levels 
of same configuration. 

In addition to the correspondence between the two sets, the program PRCBPID also calcu- 
lates the possible set of S'Lvr's for each single J-level in above group. In the set of SLns, the 
total spin is fixed while the angular momentum, L, varies. In the above example for the quintets, 
^{P,D,F,G,Hy, each J=l level is assigned to a possible set of terms, ^{P,D,F) (Table IV). 
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However, these levels can be futher identified uniquely following Hund's rule that the term with 
the larger angular momentum, L, is the lower one, i.e., the first or the lowest J—1 level should 
correspond to ^F, the second one to and the last one to ^P. 

The completeness of sets of fine structure levels with respect to the LS terms are checked. As 
mentioned above, PRCBPID determines the possible sets of S'Ltt from the target term and va- 
lence electron angular momentum of a level at its first occurance and calculates the total J- values 
of the set of LS terms. The number of these J-values, Nlv^ is compared with that of calculated 
levels, Ncal to check the completeness. For example, for the above case of ^{P, D, F,G, Hy in 
Table IV discussed above, ^P can have J — 1,2,3, can have J — 0,1,2,3,4, and so on, giving 
a total of 23 fine structure levels for this set of LS terms. The one J — level belongs to ^D, 
the three J—1 levels belong to ^{P,D,F), and so on. All 23 levels of this set are found in 
the computed levels (Table IV), thus making the computed set complete. This procedure, in 
addition to finding the link between the two diifcrcnt coupling schemes, enables an independent 
counting of the number of levels obtained, and ascertains missing or mis-identified levels. 

c) Transition probabilities 

The oscillator strengths (f-values) and transition probabilites (A-values) for bound-bound fine 
structure level transitions in Fe V are calculated for levels up to J < 8. Computations are carried 
out using STGBB of the BPRM codes. 

The f-values are initially calculated by the program STGBB with level designations given 
by Jtt only. However, the transitions may be fully described following the level identifications as 
described in the previous section. Work is in progress to identify all the transitions with proper 
quantum numbers, configurations and possible SLn^s. 

A subset of the large number of transitions has been processed with complete identifications. 
Among these transitions are those that correspond to the experimentally observed levels [23]. 
As these levels have been identified, their oscillator strenghts could be sorted out from the file of 
f-values. Another subsidiary code, PRCBPRAD, is developed to reprocess the transition proba- 
bilities where the calculated transition energies are replaced by the observed ones for improved 
accuracy. 

The computation time required for the BPRM calculations was orders of magnitude longer 
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compared to oscillator strengths calculations in LS coupling, as carried out under the OP for 
example. The time excludes that needed for creating the necessary bound state wavefunctions 
and calculating dipole matrix elements using the R-matrix package of codes. Computations are 
carried out for one or a few pairs of symmetries at a time requiring several hours of CPU time 
on the Cray T94. The memory requirement was over 30 MWords. 

4. Results and discussion 

Theoretical spectroscopic data are calculated on a large-scale with relativistic fine structure 
included in an ab initio manner, and ensuring completeness in terms of obtaining nearly all 
possible energy levels and transition probabilities for Fe V for the total angular symmetries 
considered. The results are described below. 

a) Energy Levels 

We have calculated 3,865 fine structure bound levels, with < J < 8, for Fe V. Following level 
identification, as explained in the previous section, the energy levels are arranged according to 
ascending order in energy. 

The present energies are compared with the relatively small set of experimentally observed 
levels compiled by NIST [23] in Table III. All 179 observed levels are obtained and identified. 
Asterisks attached to levels in Table III indicate an incomplete set of observed levels correspond- 
ing to the LS term. Often in experimental measurements the weak lines are not observed. The 
theoretical datasets on the other hand are usually complete. 

We find some discrepancies regarding the identification of a couple of levels in the NIST 
tabulation. The J = 2 level at 2.9395 Ry identified in the NIST table as 3d^{^P)4p{^S")2, 
from the maximum leading percentage, may have been misidentified. Present analysis for the 
completeness of a set of fine structure levels belonging to a term indicates it as an extra level for 
the given configuration and that the possible LS terms for this level are 'id^{^D2)Ap{^PDF"), 
possibly ^F°. Similarly the NIST identification for the J=3 level at 2.8968 Ry is ?,(P{^P)Ap{^D%, 
from the maximum leading percentage. Present calculations however assign the level to possible 
LS terms, 3d%'^D2)Ap{^DF''), and most likely to ^D°. 
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In the computed set of fine structure levels the observed levels are usually the ones with the 
lowest energy in each subset of Jtt. The lowest calculated levels are the 34 levels of the ground 
configuration Sd"^ of Fe V, in agreement with the observed ones. The agreement between the 
observed and calculated energies for these levels is within 1%. The calculated energies agree to 
about 1% with the measured ones for most of the observed levels. Although the energies are 
exoected to be highly accurate, but the uncertainty in the calculations is not comparable to that 
in spectroscopic observations (of the order of few wavenumbers) . 

Employing the completeness procedure the computed fine structure levels are tabulated, 
according to the two sets of cross-correlating quantum numbers: one according to the collision 
channels identified as {CfSt Lf J^TTt nl [K\s) J tt, and the other according to the complete set of 
J- values for each multiplicity (25 + 1), L and tt. A subset of the complete table of fine structure 
levels is presented in Table IV. (The complete table will be available electronically). Each set 
of levels is grouped by the possible set of LS terms followed by the levels of same configuration, 
core term, total spin multiplicity and parity, and with different J-values. The header for each 
group contains the total number of possible J-levels, Nlv, total spin multiplicity, parity, and all 
possible L values formed from the core and the outer electron. The possible J-values for each 
SLt: are given within parentheses next to each L value. 

The two sets of quantum numbers are compared. The levels that may be missing or mis- 
identified are thereby checked out. The number of computed levels, Ncal, is compared with that 
expected from angular and spin couphngs, Nlv. For most of the configurations the set of levels 
is complete except for the high lying ones. The comparison detects missing levels. An example 
is shown in the the set of 2>d?2{^D)M^{S, P, D, F, Cy in Table IV where one level with = 4 is 
missing. 

In Table IV, the effective quantum number v is specified alongwith other quantum numbers 
for each level. The consistency va. v — —i=-^ , where Et is the corresponding target energy, for 

V (E-Et) 

each set of levels may be noted. 

The possible SLns for each level are given in the last column. The levels with a single 
possible term only are uniquely defined. However, those with two or multiple term assignments 
can be defined uniquely applying Hund's rule that the higher L corresponds to the lower energy 
of same Jtt as explained in the previous section (we note that Hund's rule may not always apply 
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in cases of strong CI) . 

There are 112 levels of odd parity that we could not properly identify. Some of these levels are 
given in Table IV. These levels could be equivalent electron levels of configuration, 3p^ P")3d^ . 
The 16 LS terms of 3d\ which are '^Dl, ^P3, ^D3, ^F3, ^G3, ^H3, ^P5, ^F5, ^^5, ^D5, ^F5, 
^G5, ^J5, ^D5, ^G5, and ^55, in combination with the parent core ^P°, form 88 LS terms with 
31 singlets, 43 triplets, 13 quintets and 1 septet. The number of fine structure levels from these 
terms exceed the 112 computed levels that have not been identified. 

This new procedure of cross-correlation between two coupling schemes thus provides a pow- 
erful check on the completeness and level identification, and is expected to be of use in further 
BPRM work on complex atomic systems. 

b) Transition Probabilities 

The oscillator strengths (f-values) and transition probabilites (A-values) for fine structure level 
transitions in Fe V are obtained for J < 8. The allowed A J = 0, ±1 transitions include both the 
dipole allowed (AS* = 0, ±1) and the intercombination (AS* 7^ 0) transitions. The total number 
of computed transition probabilities is well over a million, approximately 1.46 x 10®. For most 
allowed pairs of Jtt symmetries, there are about 10^ — 10^ transitions. 

As explained in the previous section, a subset of the encoded transitions have been processed 
to present them with proper identifications. These correspond to the levels that have been 
observed. A sample of these is presented in Table V. In all of the f-values presented the calculated 
transition energy has been replaced by the observed one, using the BPRM line strengths {S) 
which are energy independent. Since measured energies in general have smaller uncertainties 
than the calculated ones, this replacement improves the accuracy of the oscillator strengths. The 
transitions among the 179 observed levels correspond to 3727 oscillator strengths. (The complete 
set of transition probabilities will be available electronically.) 

The f-values in Table V have been reordered to group the transitions of the same multiplet 
together. This enables a check on the completeness of the set of transitions. As this table 
corresponds to transitions among observed levels only, the completeness depends on the set of 
observed levels belonging to the LS terms. For the dipole allowed transitions, the LS multiplets 
are also given at the end of jj' transitions. 
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To our knowledge, no measured f- values for Fe V are available for comparison. Current NIST 
compilation [1] contains no f- values for any allowed transition. On the other hand, Fe V oscillator 
strengths for a large number of transitions were obtincd in the close coupling approximation 
under the OP [24] and the IP [25]. Both of these datasets are non-relativistic calculations 
in LS coupling and do not compute fine structure transitions. Fawcett has [26] carried out 
semi-empirical relativistic atomic structure calculations for fine structure transitions in Fe V. 
Comparison of the present /-values is made with the previous ones in Table VI, showing various 
degrees of agreement. Present values agree within 10% with those by Fawcett for a number of fine 
structure transitions of multiplets, Sd'^i^D) 3d^{^F)4p{^D°), and 3d\^D) 3d^{^P)4p{^P"), 
and the disagreement is large with other as well as with those of 3d'^{^D) — > 3(i^(^F)4p(^F°). 
The agreement of the present LS multiplets with the others is good for transitions 3d'^{^D) — > 
3d^{'^F)Ap{^F",^ D",^ P°). More detailed comparisons will be made at the completion of this 
work. 

The procedure of substitution of experimental for calculated energies provides an indication 
of uncertainties in the calculated /-values. The difference between the /-values obtained using 
the calculated transition energies and the observed ones is only a few percents (< 5%) for most of 
the allowed transitions. The difference is usually larger for the intercombination transitions which 
have lower transition probabilities. In atomic structure calculations, it is possible to rc-adjust 
eigenenergies of the Hamiltonian to match the observed ones and then use the wavefunctions to 
obtain the transition probabilities. Such a re-adjustment is not carried out in the BPRM calcu- 
lations of bound states, which are entirely ab initio, with the associated advantage of consistent 
uncertainties for most transitions considered. 

To obtain an estimate of the accuracy of the wavefunctions employed in the length and the 
velocity formulations, we plot, for example, the ^'/-values for transitions (J = ly — {J — 2)° 
and (J = 3)^ — {J — 4)° in Fig. 1. The top panel contains over 13,300 transitions between the 
pair of symmetries (J = 1)^ — {J — 2)", and the bottom panel contains over 20,200 transitions 
between the pair (J = 3)^ — ( J = 4)°. The plots show practically no dispersion for the strongest 
transitions with Q'/ ~ 5 — 10, and some dispersion around 10-20% for others with gf < 3. Up 
to gf < 0.1 the dispersion in length and velocity remains around the 10-20% level for most of 
the transitions, although the number of outlying transitions increases with decreasing gf. Given 
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the large number of points in the figures, the relatively low dispersion of gfi and gfv indicates 
that the /-values {gf divided by the statistical weight factor, 2J+1) for most of the transitions 
with gf r^l should be within 20% uncertainty. The /l's are usually more accurate than /y's 
since the asymptotic region wavefunctions are more accurately represented in the close coupling 
calculations using the R-matrix method. 

In general the intercombination transitions are weaker than the dipole allowed ones; the 
f-values can be orders of magnitude lower. The BP Hamiltonian in the present work (Eq. 2) 
does not include the two-body spin-spin and spin-other-orbit terms of Breit interaction [22]. A 
discussion of these terms is given by Mendoza etal. in a recent IP paper [27]. Their study on 
the intercombination transitions in C-like ions shows that the effect of the two-body Breit terms, 
relative to the one-body operators, decreases with Z such that for Z — 26 the computed A- values 
with and without the two-body Breit terms differ by less than 0.5 %. However, the differences 
towards the neutral end of the C-sequence is up to about 20%. It may therefore be expected that 
for Fe V the weaker intercombination f-values may also be systematically affected to a similar 
extent (the uncertainties in the dipole allowed f-values should be much less). Further studies of 
the Breit interaction in complex atoms are needed to ascertain this effect more precisely. 

Several aspects of the present work are targets for future studies, such as atomic structure 
calculations to study the effect of configuration interaction and relativistic effects on different 
types of transitions, and a detailed quantum defect analysis along interacting Rydberg series of 
levels in intermediate coupling. These studies should provide information on the accuracy of 
particular type of transitions and groups of levels, as well as address general problems in the 
analysis of complex spectra. 

5. Conclusion 

The present work is the first study of large-scale transition probabilities computed using the ac- 
curate BPRM method for a highly complex ion. Some of the results obtained herein are expected 
to form the basis for future computational spectroscopy of heretofore intractable complex atomic 
systems using efficient collision theory methods. The computational procedures developed for 
such undertakings are described, and illustrative results are presented from the ab initio Breit- 
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Pauli R-matrix calculations for Fe V. Detailed analysis for the identification of over 3,800 fine 
structure levels of Fe V is carried out using a combination of methods that include quantum de- 
fect theory. Further work on the analysis of relativistic quantum defects in intermediate coupling 
is planned. 

Following the completion of all computations and identifications, the dataset of approx- 
imately 1.5 million oscillator strengths will be described in another publication with a view 
towards astrophysical and laboratory applications. In order to complete the dataset for practical 
applications calculations are also in progress for the forbidden electric quadrupole and magnetic 
dipole transition probabilties using the atomic structure program SUPERSTRUCTURE. 

The newly acquired theoretical capability to obtain an essentially complete description of 
radiative transitions for an atomic system should enable several new advances such as: (a) the 
synthesis of highly detailed monochromatic opacity spectra [2], (b) the simulation of "quasi- 
continuum" line spectra from iron ions [28], (c) high resolution spectral diagnostics of iron in 
laboratory fusion and astrophysical sources, and (d) the analysis of experimentally measured 
spectra of complex iron ions. 
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Figure captions: 

1. Comprarison of qJl versus gfv for bound-bound fine structure level transitions in Fe V 
obtained in BPRM approximation. 
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Table I. The 19 fine strucuture levels of Fe IV in the close coupling eigenfunction expansion of 
Fe V. List of configurations, ^j, in the second sum of ^ is given below the table. 



Term 


Jt 


&(Ryd) 


Term 


Jt 


RfRvd) 


3d3(4F) 
Sd^C^F) 


3/2 


0.0 


sd^Cp) 


1/2 


0.241445 


5/2 


0.004659 


3d^CD2) 


5/2 


0.259568 


7/2 


0.010829 


3dH^D2) 
3d:^CH) 


3/2 


0.260877 


3d^{^F) 
3^3 (4p) 


9/2 


0.018231 


9/2 


0.261755 


1/2 


0.170756 


SdH^H) 


11/2 


0.266116 


3d3(4p) 


3/2 


0.172610 


Sd^i^F) 


7/2 


0.421163 


3d3(4p) 


5/2 


0.178707 


3d^^F) 
Sd^^Dl) 


5/2 


0.424684 


SdH^G) 


7/2 


0.187871 


5/2 


0.653448 


Sd^^G) 

3r/'^(2p) 


9/2 


0.194237 


Sd^^Dl) 


3/2 


0.656558 


3/2 


0.238S8S 









$j :, 3523^*^3^"^, 35^3/3^345, 3523/3^34^, 3s23/3d34d, 
3523/3^445, 3523/3^^4^, 3523/3^345, 3/3d5454d, 
3/3d54p, 3523/3d^ 3523/3^^4^, 353/3^^454^, 
3p63(i3454/, 3523p^3d^4/, 3s23/3d'^452, 353/3d^4s, 
3523/3o^^ 3p63o^^ 3/3^345240?, -is'ip^Zd'^U 



Table 11. Sample set of calculated energy levels (in -scale) of J — 2. 



E{Ry) 



'9 



E{Ry) 



-1.492946E-01 
-1.421287E-01 
-1.381112E-01 
-1.355085E-01 
-1.202324E-01 
-9.421634E-02 
-9.385297E-02 
-9.098725E-02 
-8.971136E-02 
-8.627750E-02 
-8.458038E-02 
-8.420175E-02 
-8.338890E-02 
-8.306503E-02 
-8.180410E-02 
-7.919621E-02 
-7.815480E-02 
-7.708061E-02 
-7.594121E-02 
-7.315088E-02 
-7.028267E-02 
-6.955230E-02 
-6.765945E-02 
-6.518968E-02 
-6.190890E-02 
-5.831623E-02 



2.58808 
2.65252 
2.69083 
2.71654 
2.88396 
3.25789 
3.26419 
3.31520 
3.33869 
3.40448 
3.43847 
3.44619 
3.46295 
3.46969 
3.49633 
3.55343 
3.57703 
3.60186 
3.62879 
3.69735 
3.77204 
3.79179 
3.84446 
3.91661 
4.01905 
4.14100 



-1.455232E-01 
-1.383809E-01 

-1.373195E-01 
-1.307289E-01 
-1.185430E-01 
-9.408765E-02 
-9.375868E-02 
-9.071538E-02 
-8.651911E-02 
-8.536568E-02 
-8.448129E-02 
-8.371297E-02 
-8.317358E-02 
-8.225844E-02 
-8.135593E-02 
-7.827940E-02 
-7.770277E-02 
-7.691196E-02 
-7.370981E-02 
-7.138001E-02 
-6.977748E-02 
-6.845060E-02 
-6.721368E-02 
-6.265635E-02 
-6.111942E-02 



2.62140 
2.68820 

2.69857 
2.76576 
2.90444 
3.26012 
3.26584 
3.32016 
3.39973 
3.42262 
3.44049 
3.45624 
3.46743 
3.48666 
3.50595 
3.57418 
3.58742 
3.60581 
3.68330 
3.74293 
3.78567 
3.82218 
3.85719 
3.99501 
4.04492 
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Table III. Comparison of calculated and observed energy levels of Fe V. 



Configuration 


Term 


J 


E(cal) 


E(expt) 


Q,-/4 

oa 


5 n 
u 


4 


o.o4yo 


O.oUiO 


oa 


5 n 


Q 
O 


c; cr /I o 
0.004z 


o.oUoo 


oa 


u 


Z 


c; con 


o.ouy4 


oa 


5 n 


1 


O.ODU 1 


o.oiiy 


oa 


u 


U 


O.OD^l 


o.oioz 


oa z 


3 p 


Z 


O.oz4 / 


O.Z ( ZD 


oa z 


3 p 


1 


QQQfl 

o.ov5oy 


O.ZOOD 


oa 2 


3 p 





5.3471 


5.2940 


oa 


3 F7 


D 


O.oU ( 4 




oa 


3 F7 


c 
O 


O.Olil 


a.zoOo 


Q^4 

oa 


3 TJ 

ri 


4 


O.oi4o 


OQRO 

O.zttDU 


oa z 


3 F 


4 


0.oU4o 


O.ZO i 4 


oa z 


3 TP 

r 


Q 

o 


0.oUD4 


o.zocSo 


oa z 


3 F 


Z 


O.oU ( D 




Q.r/4 

oa 




r 
O 


o.zocSi 




Q^4 




4 


O.ZDi4 


O.ZOCS4 


Q^4 

oa 


(jr 


Q 

O 


O.ZDOl 


0.z4i0 


Q j4o 

oa Z 


Cr 


4 


5.2006 


5.1798 


Q/74 
oa 


3 n 


q 


0. lyou 


o. 1 ( y4 


Oa 


3 n 


Z 


0. iy4o 


0.1 1 OZ 


Q/74 
oa 


3 n 


1 


o.iy^o 


0. 1 1 1 


oa 


1 7- 
J 


U 


0. ±c50Z 


0. ± / ±o 


oa z 




A 
U 


K 1 700 


0. ±OZU 


oa ^ 


1 £) 


2 


^ 1 '^f^'^ 

u. lOOO 


o.uy lo 


oa 


1 ^ 


Q 
O 




o.uozu 


oa 1 


3p 


2 




"4:. y'4:yo 


oa 1 


3p 


]^ 




'4:. yoyo 


3(^4 2 


3p 


Q 


4.9616 


4.9352 


3(^4 2 


3p 


4 


4 9719 


4.9460 




3p 


3 


4.9706 


A QAzlQ 




3p 


2 


4.9712 


4.9453 






4 


4.8830 


4.8636 




Id 


2 


4.6609 


4.6581 


1 







4.4302 


4.4093 


3rf^('* F)4s 


5p 


5 


3.7161 


3.7964 


3(i^('^ F)4s 


5p 


4 


3.7228 


3.8025 


3rf^('* F)4s 


5p 


3 


3.7282 


3.8077 


3rf3(4 ^)45 


5p 


2 


3.7324 


3.8116 


3rf^{'^ F)4s 


5p 


1 


3.7352 


3.8143 


Sd^C^ F)4s 


3p 


4 


3.6222 


3.7194 


3d^{'^ F)4s 


3p 


3 


3.6311 


3.7277 


3d^{'^ F)4s 


3p 


2 


3.6381 


3.7344 


3rf^{'^P)4s 


5p 


3 


3.5483 


3.6402 


3rf3^P)4s 


5p 


2 


3.5532 


3.6453 


3rf^(4p)45 


5p 


1 


3.5554 


3.6475 


3rf^(2G)45 


3G 


5 


3.5090 


3.6038 


3^^(2^)45 


3G 


4 


3.5130 


3.6076 


3d^{'^G)4s 


3G 


3 


3.5158 


3.6101 


3rf3{'*P)4s 


3p 


2 


3.4595 


3.5663 


3rf^^P)45 


3p 


1 


3.4677 


3.5738 


3rf^^P)45 


3p 





3.4706 


3.5763 


3d;^(^G)4s 


iG 


4 


3.4636 


3.5673 


3d^(^P)4s 


3p 


2* 


3.4330 


3.5583 


3d3(2p)4s 


3p 


1* 


3.4515 


3.5575 


3<i3(2_D2)4s 


3d 


3 


3.4337 


3.5399 


3^3(21)2)43 


3d 


2 


3.4528 


3.5394 


3d^{'^D2)4s 


3d 


1 


3.4405 


3.5468 


3d^{^H)4s 


3i? 


6 


3.4364 


3.5346 


3d^{^H)4s 


3if 


5 


3.4392 


3.5370 
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Table III. continues. 



Configuration Term 


J 


E(cal) 


E(cxpt) 






3 TT 


4 


3.4398 


3.5377 


3a^ 




1 D 

^P 


1 


3.4015 


3.5131 


o J'^ 


■^D2)4s 




2 


3.3877 


3.5027 


o J'^ 
3a^ 


il)4s 




5 


3.3894 


3.4965 




^F)4s 


3 j-1 


4 


3.2711 


3.3841 




7-l\ 4 

^F)4s 


3 7-1 


3 


3.2695 


3.3822 




7-l\ 4 

''FjAs 


3 7-1 

^F 


2 


3.2682 


3.3806 


3a 




1 7-1 
P 


3 


3.2280 


3.3468 


3a^ 


'*F)4p 


5qo 


6 


3.1534 


3.1594 


3d^ 


(■*F)4p 


5qo 


5 


3.1631 


3.1700 


3d^ 


(■'F)4p 


5qo 


4 


3.1710 


3.1787 




■*F)4p 


5qo 


3 


3.1773 


3.1858 


3d^ 


■*F)4p 


5qo 


2 


3.1821 


3.1912 


3d^ 


(■'F)4p 


5 j^o 


1* 


3.1495 


3.1644 


3(i 


(*'F)4p 


5£)0 


4 


3.1359 


3.1498 


3a^ 


(*'F)4p 


5£)0 


3 


3.1433 


3.1559 


3d^ 


(*'F)4p 


5£)0 


2 


3.1479 


3.1609 


3a 


(^F)4p 




1 


3.1408 


3.1540 


O 7'-i 

3a'' 


'*F)4p 


''D° 





3.1469 


3.1565 


O 7'-i 

3(i 


'•^Dl)4s 


3 7-1 


3 


3.0074 


3.1581 


O 7'^ 

3rf 


Dl)4s 


3 7-1 


2 


3.0058 


3.1564 


O 7'^ 




3 7-1 

''D 


1 


3.0047 


3.1551 


3d^ 


(■'F)4p 




5* 


3.1191 


3.1343 


3d^ 


(■'F)4p 




4* 


3.1243 


3.1391 


O 7'^ 


(*'F)4p 






3.1.300 


3.1443 


O 7'^ 

3(i 


(*'F)4p 


;> 


2* 


3.1350 


3.1496 


O 7'^ 

3d' 


(^F)4p 


'^D° 


3 


3.1170 


3.1331 


O 7'^ 


(*'F)4p 




2 


3.1258 


3.1401 


O j3 

3(i'^ 


*F)4p 


3£)0 


1 


3.1309 


3.1439 


3a'' 


^_Dl)4s 


1 7-1 


2 


2.9636 


3.1210 


3a^ 


(■*F)4p 




5 


3.0737 


3.0973 


3a^ 


(■*F)4p 




4 


3.0812 


3.1035 


O j3 

3(i^ 


■*F)4p 


3 


3 


3.0872 


3.1083 


o J'^ 
3d^ 


■*F)4p 


3 7-10 

o pw 


4 


3.0449 


3.0716 


O j'^ 

3d 


(■'F)4p 


3 JT^O 


3 


3.0520 


3.0779 


3(1 


(''F)4p 


5^o 


2 


3.0586 


3.0836 


o j3 

Sd'' 


■*P)4p 




3 


2.9842 


3.0078 


o j3 

3(i'' 


^P)Ap 


5^o 


2 


2.9903 


3.0151 


o j3 

3(i'' 




po 


1 


2.9937 


3.0195 


o 73 


*P)ip 


''D° 


4 


2.9526 


2.9792 


o j3 

3d^ 


(■*P)4p 




3 


2.9600 


2.9883 


oa 


j4p 


LI 


2 


2.9733 


2.9912 


3d^ 


('lp)4p 




1 


2.9747 


3.0058 


3d^ 


('lp)4p 


5D° 





2.9797 


3.0094 


3d^ 


(■lp)4p 


Z po 


2 


2.9622 


3.0038 


3d^ 


(■lp)4p 


3po 


1 


2.9635 


2.9911 


3d^ 


(■*P)4p 


3 po 





2.9676 


2.9941 


3d^ 


(2G)4p 


■Sfjo 


6 


2.9484 


2.9739 


3d^ 


(2G)4p 




5 


2.9591 


2.9863 


3d'' 


(2G)4p 


i^o 


4 


2.9655 


2.9942 


3d^ 


(2G)4p 




5 


2.9289 


2.9613 


3d^ 


(2G)4p 


■Aqo 


4 


2.9353 


2.9662 


3d^ 


(2G)4p 




3 


2.9412 


2.9726 


3d^ 


(2G)4p 


3 po 


4 


2.9238 


2.9583 


3d^ 


(2G)4p 


3 po 


3 


2.9262 


2.9540 


3d* 


(2G)4p 


3 po 


2 


2.9329 


2.9567 


3^3 


(2p)4p 


3po 


1 


2.9326 


2.9439 


3d3 


(2p)4p 


3po 





2.9458 


2.9413 


3d3 


[2G)4p 


IQO 


4 


2.9068 


2.9430 



Table III. continues. 



Configuration 


Term 


J 


E(cal) 


E(expt) 


3dr 




5 go 1 


2 


2.9231 


2.9395 






1 po 


3 


2.9067 


2.9382 






•'S° 


2 


2.8995 


2.9379 






3 po 


2 


2.8605 


2.8668 






Ijjo 


5 


2.9006 


2.9354 




P)4p 


3p,o 


3 


2.8895 


2.9117 




^P)4p 


3p)0 


2 


2.8944 


2.9169 


o j3 


^P)4p 


3p,o 


1 


2.8790 


2.9274 


o j3 

3d^ 




3jjo 


6 


2.8793 


2.9143 






3jjo 


5 


2.8830 


2.9180 




(0 TT\ A 


3fjo 


4 


2.8834 


2.9189 




■^D2)4p 


1 po 


1 


2.9031 


2.9073 






3po 


4* 


2.8653 


2.8922 




^D2)4p 


3 po 


2* 


2.9144 


2.9055 


O 73 


■^P)4p 


3 go 


1 


2.8621 


2.9052 




('*P)4p 


Sp)0 


3 


2.8593 


2.9030 


O 73 


*P)4p 


3p)0 


2 


2.8749 


2.8991 


o j3 




3p)o 


1 


2.8478 


2.8991 


O 73 

3d-^ 


^9 7-l\ < 

(^P)4p 


3p)0 1 


3* 


2.8718 


2.8968 


o j3 


y'^ H)Ap 


3 JO 


7 


2.8414 


2.8780 


o i3 


('? TT\ A 

^H)4p 


3 JO 


6 


2.8487 


2.8872 


o (3 


(^H)Ap 


3 JO 


5 


2.8537 


2.8938 


O 73 


(^D2)Ap 


3j5o 


3 


2.8366 


2.8713 


o j3 

3 a' 


(^D2)4p 


3j5o 


2 


2.8325 


2.8761 


3(i'^ 


('^D2)4p 


3j5o 


1 


2.8842 


2.8826 


O 73 


(^H)'ip 


iG° 


4 


2.8325 


2.8746 




(^H)ip 


Ijjo 


5 


2.8244 


2.8696 


O 73 

3d 


(^Di)Ap 


3po 


1* 


2.8313 


2.8652 


O 73 


^D2)4p 


3 po 


0* 


2.8300 


2.8623 


o i3 

3d-^ 


^D2)4p 


1 j;io 


3 


2.8238 


2.8593 




{1 TT\ A 

''H)Ap 


3 go 


5 


2.8039 


2.8496 


o j3 


<9 TT\ A 

^H)ip 


3 go 


4 


2.8035 


2.8483 


o j3 


^9 7T\ A 

■'Hjip 


3qo 


3 


2.8049 


2.8476 


o j3 


^9 7-7-\ A 

''HjAp 


Ijo 


6 


2.8012 


2.8489 


n j3 


^P)ip 


3 go 


1 


2.7705 


2.8282 


o j3 


''^D2)Ap 


Ijpo 


2 


2.7815 


2.8184 


O 73 


11 T~t\ A 


1 JJO 


1 


2.7822 


2.8161 


o j3 


^9 7-l\ A 

^P)4p 


3 po 


4 


2.7083 


2.7556 




^P)4p 


3 po 


3 


2.7105 


2.7577 


n j3 


^9 7-l\ A 

'^P)4p 


3 po 


2 


2.7109 


2.7585 


O 7 -t 


/9 T-^^ J 

;^p)4p 


'-''G° 


5 


2.6652 


2.7150 




(^P)4p 


3qo 


4 


2.6672 


2.7190 


o j3 


■^P)4p 


3qo 


3 


2.6695 


2.7229 


O 73 


(2p)4p 


3p,o 


3 


2.6648 


2.7129 


Q^3 


_r j4p 


3p,o 


2 


2.6705 


2.7050 




(2p)4p 


''•D° 


1 


2.6559 


2.7003 




(2p)4p 


Ijjo 


2 


2.6592 


2.7097 


3d-^ 


(2p)4p 


iG° 


4 


2.6165 


2.6775 


3d-^ 


(2p)4p 


1 jT^o 


3 


2.6172 


2.6742 


3d-^ 


(2Dl)4p 


3j)o 


3 


2.4918 


2.5249 


3d-^ 


(2Dl)4p 


3j)o 


2 


2.4938 


2.5278 


3d-^ 


(2Dl)4p 


3j)o 


1 


2.4936 


2.5285 


3d-^ 


(2Dl)4p 


lj)0 


2 


2.4581 


2.5074 


3d-^ 


(2Dl)4p 


3 po 


4 


2.4499 


2.4876 


3d-^ 


(2Dl)4p 


3 J?o 


3 


2.4532 


2.4935 


3d-^ 


(2Dl)4p 


3 po 


2 


2.4510 


2.4938 


3d^ 


(^Dljlp 


3 po 


2 


2.4196 


2.4649 


3d3 


(2Dl)4p 


3po 


1 


2.4141 


2.4580 


3d3 


(2Dl)4p 


3po 





2.4116 


2.4546 


3d3 


(2Dl)4p 


1 po 


3 


2.4058 


2.4518 


3(i3 


(2Dl)4p 


1 po 


1 


2.3394 


2.3924 



Table IV. Sample table of calculated and identified fine strucuture energy levels of Fe V. Nlv=total 
number of levels expected for the possible LS terms (specified as 2S + 1, tt, and set of L with 
J -values within paratheses), formed from the target term and I of the outer electron, and Ncal — 
number of calculated levels. SLir in last column— possible LS terms for each level. 



Ct 




StLtTTt Jt 


nl 


J 


E(cal) 




SL-K 


l\ IV — 




5,e: F ( 5 4 3 2 1 ) 














3d3 






'4Fe' 


1 3/2 


4s 


1 


-3.73515EH-00 


2.59 


5 


F e 


3d3 






4Fe 


1 5/2 


4s 


2 


-3.73238E+00 


2.59 


5 


F e 


oao 






4Fe 


1 5/2 


Ac- 
4S 


o 
O 




z.oy 


5 


F e 


3d3 






4Fe 


1 7/2 


4s 


4 


-3.72275E+00 


2.59 


5 


F e 


3d3 






;4Fc 


1 9/2 


4s 


5 


-3.71610E+00 


2.59 


5 


F e 


Ncal= 


5 


: set complete 














Nlv= 


3, 


3,e: F ( 4 3 2 ) 














3d3 




1 




1 3/2 


4s 


2 


-3.63808E+00 


2.62 


3 


F e 


3d3 






:Si 


1 7/2 


4s 


3 


-3.63107E+00 


2.62 


3 


F e 


3d3 




1 


;4Fc: 


1 9/2 


4s 


4 


-3.62225E+00 


2.62 


3 


F e 


Ncal= 


3 


: set complete 














Nlv= 


3, 


l,o: P ( 1 


)D(2) 


F(3; 


1 










3d3 1 




( 


2De 


) 3/2 


4p 


2 


-2.45812E+00 


2.83 


1 


D o 


3d3 1 




( 


'2De' 


) 5/2 


4p 


3 


-2.40581E+00 


2.86 


1 


F o 


3d3 1 






'2De 


) 5/2 


4p 


1 


-2.33944E+00 


2.89 


1 


P o 


Ncal= 


3 


: set complete 















Nlv= 23, 5,e: P(321)D(43210)F(54321)G(65432)H(7654 3) 



3d3 1 


'4Fe^ 


1 3/2 


4d 


3 


-2.37021E+00 


3.25 


5 


PDFGH e 


3d3 ( 


:4Fe; 


1 5/2 


4d 


4 


-2.36647E+00 


3.25 


5 


DFGH e 


3d3 1 


'4Fc^ 


1 5/2 


4d 


5 


-2.36189E+00 


3.25 


5 


FGH e 


3d3 1 


:4Fc; 


1 3/2 


4d 


1 


-2.35988E+00 


3.25 


5 


PDF e 


3d3 1 


'4Fc' 


1 7/2 


4d 


6 


-2.35651E+00 


3.25 


5 


GH e 


3d3 1 


*4Fc 


1 5/2 
1 9/2 


4d 


2 


-2.35541E+00 


3.26 


5 


PDFG e 


3d3 1 


*4Fc 


4d 


7 


-2.35041E+00 


3.25 


5 


H e 


3d3 1 


*4Fe* 


1 5/2 


4d 


1 


-2.34932E+00 


3.26 


5 


PDF e 


3d3 1 


*4FcJ 


1 9/2 


4d 


3 


-2.34736E+00 


3.25 


5 


PDFGH e 


3d3 1 


'4Fe' 


1 7/2 


4d 


2 


-2.34633E+00 


3.26 


5 


PDFG c 


3d3 1 


*4Fe^ 


1 3/2 


4d 


2 


-2.34397E+00 


3.27 


5 


PDFG c 


3d3 1 


*4Fc 


1 7/2 


4d 


4 


-2.34329E+00 


3.26 


5 


DFGH c 


3d3 1 


*4Fc 


1 5/2 


4d 


3 


-2.34092E+00 


3.26 


5 


PDFGH c 


3d3 ( 


*4FcJ 


1 9/2 


4d 


3 


-2.33989E+00 


3.26 


5 


PDFGH c 


3d3 1 


'4Fe' 


1 9/2 


4d 


5 


-2.33822E+00 


3.26 


5 


FGH c 


3d3 ( 


*4Fe 


1 7/2 


4d 


5 


-2.33234E+00 


3.27 


5 


FGH e 


3d3 ( 


*4Fe 


1 9/2 


4d 


6 


-2.32699E+00 


3.26 


5 


GH e 


3d3 ( 


;4Fe' 


1 3/2 


4d 


4 


-2.28772E+00 


3.31 


5 


DFGH e 


3d3 1 


'4Fe' 


1 3/2 


4d 





-2.28673E+00 


3.31 


5 


D e 


3d3 ( 


*4Fe 


1 7/2 


4d 


1 


-2.28265E+00 


3.30 


5 


PDF e 


3d3 ( 


*4Fe 


1 9/2 


4d 


2 


-2.27468E+00 


3.30 


5 


PDFG e 


3d3 ( 


*4Fe 


1 7/2 


4d 


3 


-2.26346E+00 


3.31 


5 


PDFGH e 


3d3 ( 


*4Fc 


) 9/2 


4(1 


4 


-2.25835E+00 


3.31 


5 


DFGH c 



Ncal= 2:i : set (-(jmplcic 



Nlv= 15, 3,e: P(210)D(321)F(432)G(543)H(6 5 4) 



3d3 ( 


'4Fe' 


1 3/2 


4d 


1 


-2.35634E+00 


3.26 


3 


PD e 


3d3 ( 


*4Fe 


1 5/2 


4d 


2 


-2.35219E+00 


3.25 


3 


PDF e 


3d3 1 


*4Fe 


1 7/2 


4d 


3 


-2.34872E+00 


3.25 


3 


DFG 6 


3d3 1 


*4Fe; 


1 5/2 


4d 


4 


-2.33701E+00 


3.27 


3 


FGH e 


3d3 1 


'4Fe' 


1 5/2 


4d 


5 


-2.28113E+00 


3.31 


3 


GH e 


3d3 1 


*4Fe 


1 3/2 


4d 


3 


-2.28060E+00 


3.31 


3 


DFG e 


3d3 1 


*4Fe 


1 7/2 


4d 


4 


-2.27417E+00 


3.31 


3 


FGH e 


3d3 1 


*4Fe 


1 9/2 


4d 


6 


-2.27323E+00 


3.30 


3 


H 6 


3d3 1 


*4Fc; 


1 3/2 


4d 


2 


-2.26789E+00 


3.32 


3 


PDF e 


3d3 1 


'4Fe' 


1 7/2 


4d 


5 


-2.26632E+00 


3.31 


3 


GH e 


3d3 1 


'4Fe 


1 5/2 


4d 





-2.24585E+00 


3.33 


3 


P 6 


3d3 1 


*4Fe; 


1 5/2 


4d 


1 


-2.24483E+00 


3.33 


3 


PD e 


3d3 1 


'4Fc' 


1 7/2 


4d 


2 


-2.24278E+00 


3.33 


3 


PDF e 


3d3 1 


*4Fc 


1 7/2 


4d 


3 


-2.23973E+00 


3.33 


3 


DFG e 


3d3 ( 


*4Fe; 


1 9/2 


4d 


4 


-2.23571E+00 


3.33 


3 


FGH e 



Ncal= 15 : set complete 



Nlv= 13, 5,e: P(321)D(43210)F(54321) 








3d3 ( 


'4Pc^ 


I 1/2 


4d 


1 


-2.16405E+00 


3.27 


5 


PDF c 


3d3 ( 


'4Pc' 


1 3/2 


4d 


2 


-2.16298E+00 


3.27 


5 


PDF c 


3d3 ( 


'4Pe' 


) 3/2 


4d 


3 


-2.16165E+00 


3.27 


5 


PDF e 


3d3 ( 


'4Pe* 


) 3/2 


4d 


1 


-2.15901E+00 


3.27 


5 


PDF e 


3d3 ( 


'4Pe* 


) 3/2 


4d 


4 


-2.15825E+00 


3.27 


5 


DF e 


3d3 ( 


'4Pe* 


) 5/2 


4d 


2 


-2.15694E+00 


3.27 


5 


PDF e 


3d3 ( 


'4Pe* 


I 5/2 


4d 


3 


-2.15521E+00 


3.27 


5 


PDF e 


3d3 ( 


'4Pe* 


1 5/2 


4d 


5 


-2.15477E+00 




5 


F e 



Table IV. continues. 



Ct 



StLtTTt 



Jt 



nl 



E(cal) 



SLn 



Nlv= 13, 5,e: P ( 3 2 1 

3d3 (4Pe) 

3d3 (4Pe) 

3d3 (4Pc) 

3d3 (4Pc) 

3d3 (4Pe) 

Ncal= 13 : set complete 



)D(43210)F(54321) 

1/2 4d 3 -2.13511E+00 3.29 

5/2 4d 1 -2.13294E+00 3.29 

1/2 4d 2 -2.11203E+00 3.31 

5/2 4d 4 -2.09728E+00 3.31 

3/2 4d -2.08337E+00 3.33 



PDF e 
PDF e 
PDF e 
DF e 
D e 



Nlv= 

3d3 
3d3 
3d3 
3d3 
3d3 
3d3 
3d3 
3d3 
3d3 
Ncal= 



9, 3,e: P(210)D(3 21)F{432) 



■4Pc" 


1 5/2 


4d 


2 


-2.13414E+00 


3.29 


3 


PDF e 


■4Pc' 


) 3/2 


4d 


3 


-2.10788E+00 


3.31 


3 


DF e 


'4Pe' 


) 3/2 


4d 


1 


-2.08179E+00 


3.33 


3 


PD e 


'4Pe' 


) 5/2 


4d 


2 


-2.07934E+00 


3.46 


3 


PDF e 


'4Pe' 


) 5/2 


4d 


1 


-2.07801E+00 


3.33 


3 


PD e 


'4Pc* 


) 5/2 


4d 


3 


-2.07780E+00 


3.33 


3 


DF e 


'4Pc* 


1 5/2 


4d 





-2.07641E+00 


3.33 


3 


P e 


'4Pe* 


I 5/2 


4d 


4 


-2.07512E+00 


3.33 


3 


F e 


'4Pe* 


1 3/2 


4d 


2 


-2.04510E+00 


3.36 


3 


PDF e 



9 : set complete 



Nlv= 15, 3,e: D{321)F{432)G(543)H(654)I(765) 



3d3 ( 


'2Gc; 


1 7/2 


4d 


3 


-2.16065E+00 


3.26 


3 


DFG e 


3d3 ( 


;2Ge; 


1 7/2 


4d 


5 


-2.15935E+00 


3.26 


3 


GHI e 


3d3 ( 


'2Ge; 


1 7/2 


4d 


6 


-2.15443E+00 


3.27 


3 


HI e 


3d3 ( 


'2Ge^ 


1 7/2 


4d 


4 


-2.15035E+00 


3.27 


3 


FGH e 


3d3 ( 


'2Ge^ 


1 9/2 


4d 


7 


-2.14872E+00 


3.27 


3 


I e 


3d3 ( 


'2Ge* 


1 9/2 


4d 


3 


-2.14865E+00 


3.27 


3 


DFG e 


3d3 ( 


'2Gc' 


1 9/2 


4d 


5 


-2.14663E+00 


3.27 


3 


GHI e 


3d3 ( 


'2Gc* 


1 9/2 


4d 


5 


-2.13490E+00 


3.28 


3 


GHI e 


3d3 ( 


'2Gc* 


1 9/2 


4d 


6 


-2.13050E+00 


3.30 


3 


HI e 


3d3 ( 


'2Gc; 


1 9/2 


4d 


2 


-2.11451E+00 


3.29 


3 


DF e 


3d3 ( 


;2Ge; 


1 7/2 


4d 


1 


-2.10743E+00 


3.30 


3 


D e 


3d3 ( 


'2Ge' 


1 7/2 


4d 


2 


-2.10505E+00 


3.30 


3 


DF e 


3d3 ( 


'2Ge 


1 7/2 


4d 


3 


-2.09600E+00 


3.31 


3 


DFG e 


3d3 1 


'2Ge 


1 9/2 


4d 


4 


-2.08773E+00 


3.31 


3 


FGH e 


3d3 ( 


*2Ge 


1 9/2 


4d 


4 


-1.95461E+00 


3.41 


3 


FGH e 



Ncal= 15 : set complete 



Nlv= 5, l,e: 



D(2)F(3)G{4)H(5)I(6) 



3d3 


(2Ge' 


1 7/2 


4d 


5 


-2.15108E+00 


3.27 


1 


H c 


3d3 


(2Ge^ 


1 7/2 


4d 


4 


-2.13926E+00 


3.27 


1 


G c 


3d3 


(2Gc* 


1 9/2 


4d 


3 


-2.11077E+00 


3.29 


1 


F e 


3d3 


(2Gc* 


1 7/2 


4d 


6 


-2.10182E+00 


3.27 


1 


I e 


3d3 


(2Ge* 


1 7/2 


4d 


2 


-2.09282E+00 


3.30 


1 


D e 



Ncal= 5 : set complete 



Nlv= 9, 3,e: 

3d3 
3d3 
3d3 
3d3 
3d3 
3d3 
3d3 
3d3 
3d3 

Ncal= 9 : set complete 




D(321)F(432) 



-2.11449E+00 
-2.10214E+00 

-2.08472E+00 
-2.08116E+00 
-2.07662E+00 
-2.05165E+00 
-1.94013E+00 
-1.93293E+00 
-1.89853E+00 



3.26 
3.27 

3.28 
3.28 
3.29 
3.30 
3.39 
3.39 
3.42 



PD e 

F e 
PDF e 
DF e 
PDF e 
DF e 
P e 
PD e 
PDF e 



Nlv= 3, l,e: P(1)D(2)F(3) 

3d3 (2Pe) 3/2 4d 1 -2.08815E+00 3.26 1 P e 

3d3 (2Pe) 1/2 4d 2 -1.95699E+00 3.41 1 D e 

3d3 (2Pe) 3/2 4d 3 -1.94667E+00 3.38 1 F e 

Ncal= 3 : set complete 



Nlv= 5, 1,' 



S(0)P{1)D(2)F(3)G(4) 



3d3 2 ( 


'2Dc) 


3/2 


4d 


3 


-2.06726E+00 


3.28 


1 


F c 


3d3 2 ( 


'2De) 


3/2 


4d 
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1 


P e 


3d3 2 1 


'2De) 


5/2 


4d 
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-1.92280E+00 


3.38 


1 


D e 


3d3 2 ( 


*2De) 


3/2 


4d 
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-1.88850E+00 


3.46 


1 


G e 


3d3 2 ( 


'2Dc) 


5/2 


4d 





-1.84426E+00 


3.45 


1 


S e 



Ncal= 5 : set complete 
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Table IV. continues. 



Ct 



StLtTTt 



nl 



E(cal) 



Nlv= 13, 

3d3 2 
3d3 2 
3d3 2 
3d3 2 
3d3 2 
3d3 2 
3d3 2 
3d3 2 
3d3 2 
3d3 2 
3d3 2 
3d3 2 

Ncal= 12 



S(1)P(210)D(321)F(4 3 2)G(54 3) 




-1.05446E+00 4.36 

-1.04247E+00 4.38 

-1.04179E+00 4.38 

-1.03332E+00 4.40 

-1.03014E+00 4.40 

-1.02889E+00 4.40 

-1.02130E+00 4.42 

-1.01682E+00 4.43 

-1.01420E+00 4.43 

-1.00936E+00 4.44 

-9.92578E-01 4.47 

-9.81500E-01 4.49 



3 SPD e 
3 G e 
3 DFG e 
3 SPD e 
3 FG e 
3 PDF e 
3 P c 
3 SPD e 
3 DFG e 
3 PDF e 
3 PDF e 
3 DFG e 



level missing: 4 



Odd Parity, Equivalent-electron & unidentifiable levels (see discussion) 



3p53d5 


1 


-1.23559E+00 


3p53d5 
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-1.22406E+00 


3p53d5 
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-1.20599E-I-00 


3p53d5 


3 


-1.20156E+00 


3p53d5 


2 
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3p53d5 


5 


-1.12662E+00 
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-1.12076E+00 
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-1.11088E+00 


3p53d5 
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-1.06928E-I-00 
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-1.02720E+00 
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-9.99298E-01 


3p53d5 


1 


-9.99223E-01 


3p53d5 
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-9.98823E-01 


3p53d5 


5 


-9.92970E-01 
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Table V. Transition probabilties of Fe V among observed fine structure levels. 
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Table V. continues. 
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Table VI. Comparison of present f -values with the earlier ones. 
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a Fawcett (1989), b Butler (TOPbase), c Bautista (1996) 
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